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Gradual Damping of Solitary Waves 

By Garbis H. Keulegan 



This paper treats the problem of the damping by viscous action of translation waves. 
A short exposition is given of Boussinesq's boundary layer theory for wave motion, and 
expressions for the damping of rectangular and solitary waves are derived. Scott Russell's 
experimental results for solitary waves are compared with the theory, and satisfactory 
agreement is found to exist. This fact makes it legitimate to apply the formulas developed 
to correct in model tests on harbors or in other tests of a like nature for the dissipative effects 
that occur in shallow-water waves. 



I. List of Symbols 



A = 
B= 

e = 

dE l /dt = 

dE 2 /dt= 

8E X = 

F( ) = 

0= 

h or hi = 

ho\ = 

H-- 

k~- 

K-- 

, L- 

AT: 
P= 

r- 



z-- 
w, v f w= 

u 0y ^o, w -- 



X- 
Y- 
Z- 



■ohlH. 

breadth of channel. 
: base of natural logarithms, 
rate of dissipation of energy of the wave. 
: rate of decrease of energy of the wave. 
: energy in column of unit length and of width dx. 
; function of. 
-acceleration of gravity. 
= height of wave above undisturbed surface, 
initial height of wave above undisturbed surface, 
-depth of undisturbed water. 
= constant defined by eq 35. 
: constant defined by eq 44. 
= length of wave. 
-kx. 

= an integral denned by eq 41. 
= the pressure in the liquid. 
= a variable. 
-a variable, also the distance travelled by the 

wave. 
: time. 
= coordinate parallel to solid surface in direction 

of motion. 
: coordinate parallel to solid surface at right angles 

to direction of motion. 
^ coordinate perpendicular to solid surface. 
= velocity components of particles in the boundary 

layer parallel to x, y, z, respectively. 
= velocity components of particles in the potential 

region parallel to x, y, z, respectively. 
-- particle velocity in the potential region, 
^-component of body force per unit mass due to 

gravity. 
= ?/-component of body force per unit mass due to 

gravity. 
= z-component of body force per unit mass due to 

gravity. 



Y( ) = function of. 

a = a variable, also a damping coefficient. 

/3 = a variable. 

/c = a constant defined by eq 17. 

M = viscosity of the water. 

v = kinematic viscosity of the water. 

p = density of the water. 

<f> = & dissipation function. 

co = velocity of wave propagation. 

II. Introduction 

In the propagation of translation waves in still 
water, the main source of dissipation of the energy 
is the thin layer of liquid next to the solid bound- 
aries of the channel; that is, in the boundary layer. 
The motion in the layer is laminar, and the de- 
termination of the velocities in it can be made by 
resorting to the ordinary equations of viscous 
motion. The motion of the liquid outside of the 
layer is potential, and there the values of the 
velocities are exactly those that would be obtained 
from the ordinary theory, assuming that the lami- 
nar layer is absent. The errors involved in the 
use of this method can be ignored, since the varia- 
tion of pressure in a cross section of the laminar 
layer is negligible. In fact, the gradual decrease of 
height of a solitary wave, the only wave of transla- 
tion, which under ideal conditions, travels without 
change of form, can be explained on the basis of 
this dissipation, at least, when the wave is of 
moderate proportions. The purpose of this paper 
is to develop a method of computing this energy 
dissipation and to compare the results with ex- 



Gradual Damping of Solitary Waves 



487 



perimental data as a means of verifying the 
accuracy of the theory. 

III. Dynamical Equation of the Laminar 
Layer 

The theory of the laminar layers with reference 
to wave motion was initiated by Boussinesq, and 
the essential elements of his method will now be 
given [I]. 1 

Consider the unsteady motion of a liquid near 
a flat solid boundary. For the sake of generality 
it vvill be supposed that the motion is two-dimen- 
sional. We select the rectangular system of axes 
x 7 y, z, where z is normal to the solid surface. The 
velocity components of particles along these axes 
are u, v, and 0, respectively. At the solid bound- 
ary, that is, at z=0, u and v vanish. In the motion 
of waves u, du/dx, and similar quantities are small, 
and hence their products may be neglected. 
Owing to the thinness of the layers and their great 
lengths, only shears parallel to the solid boundaries 
are of consequence. Thus, the equation of motion 
for the liquid in the layer may be written as 



du 
d~t~ 

dv 
dt' 



p dx di 2 

—pdy +Y+v ^ 



ds 2 



(1) 



and 



= - 



1 dp 



p dz 



;+Z, 



where v is the kinematic viscosity, p is the density 
of the liquid, p is the pressure, and X, Y, Z are the 
components of the body forces per unit mass due 
to gravity. The third equation indicates that in 
the boundary layer the pressure p, the part due 
to disturbances, does not vary along the normals 
to the boundary. The pressure, p, in the viscous 
layer is the same function of x and y as in the 
potential region next to the viscous layer. 

Denoting the limiting values of u and v at the 
boundary of the layer by u and % these last being 
the velocities in the potential region, we have for 
this region 



du ^ 1 dp y \ 
dt " pdx 1- '' 

dt pdy^ 1 ' 



(2) 



Subtracting these equations from the corre- 
sponding ones in eql, we obtain 



d(u— Uq) _ d 2 (u—u ) 



dt 



dz 2 



and 



d{v—v Q ) d 2 (v—v ) 
~dT =v ~dz^' 

with the boundary conditions 

u—u =—u z=Q 

U — U = 2= oo 



and 



v—v Q =—v z=0 

V — V =0 Z= oo 



(3) 



(4) 



(5) 



(6) 



These are the basic boundary equations as de- 
rived by Boussinesq. In the form given above 
they are applicable to a progressive wave of 
oscillation in a channel or to a standing wave in a 
rectangular basin. 

The limiting value, z= &> } appearing in the 
boundary conditions is introduced for convenience 
in the mathematical analysis. Physically the 
boundary layer is very thin, and the velocities 
u and v attain the values u and v , respectively, 
at small finite distances away from the solid 
boundaries. Moreover, if the waves examined 
are progressive oscillatory waves moving in a 
rectangular channel and the layer next to the 
bottom is considered, the axis of z is taken vertical 
and the x axis is in the same direction as the 
channel axis. For this case z?=0, and u is a 
function of x, z, and t. When the boundary layer 
of one of the vertical walls is considered, the axis 
of z is taken normal to the wall, the axis of y is 
vertical, and, as before, the axis of x is parallel to 
the channel axis. Now u and v are the velocities 
in the potential region outside of the layer and are 
functions of x, y, and t only. The case of the 
standing waves in a rectangular channel or in a 
rectangular basin is interpreted in the same 
manner. 

As we shall be interested primarily in transla- 
tion waves, we need consider only 



1 Figures in brackets indicate the literature references at the end of this 
paper. 



d(u—u ) _ d 2 (u—Up) 
~dT ~ v dz 2 ' 



(7) 



488 



Journal of Research 



with the boundary conditions 

u— u = — u 2=0 

u—u = 2 = oo 



(8) 



Now u is the uniform velocity in the cross section 
and is a function of x and t, assuming that the 
wave travels in the direction of x positive or x 
negative. 

In the approximate theory of waves, if h is the 
superelevation of the wave surface, measured from 
the undisturbed water surface, and H is the depth 
of the water in the channel, for a wave traveling in 
the direction of x positive: 



u -- 



h 



(9) 



where co is the velocity of wave propagation given 
by 

co=V<7#. (10) 

For a wave of translation having any arbitrary 
form but a small wave height, and for a solitary 
wave, independently of height, for waves traveling 
in the direction of x positive, we have 



and 



h=h(x — oot), 
u=u (x—ut). 



(ID 
(12) 



The meaning of these expressions is that the 
shape of waves does not change with translation. 
Obviously there are two methods available for the 
description of the waves. We may observe the 
variation of the wave height with time at a fixed 
point x(x=0), or we may observe the variation of 
the wave height with x at a fixed time t(t=0). 

We select the second method of observation, and 
then the main problem resolves itself into de- 
termining the velocities in the boundary layer 
in the entire length of the wave at a given time, t. 
From this point of view, the basic equation now is, 
replacing d/d£ in eq 7 by — a>(d/dx) 



d(u—u __ v b 2 (u—u ) 



dx 



d2 2 



ith the boundary conditions 

u— u = — u Q 2 = 

U — U = 2= co 
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(13) 



(14) 



In general, for t constant, it will be 


assumes 


thai 


u =F(x). 




(15) 


IV. Distribution of the Velocities in the Vis- 
cous Layer. 


The expression 






u—u =—e cos K{a—x)- 


-V§} 


(16) 


where k and a are independent of 2 and x, as is 
well known, is a solution of eq 13 and satisfies the 
second equation of the boundary conditions, eq 14. 
Putting 2=0, and remembering that u=0 at 
2=0, we obtain 


u =cos [k(qi — x)]. 




(17) 


Thus we conclude that when the velocity of u 
in the potential region of the water is periodic and 
is given by eq 17, velocity u in the laminar layer 



is given by eq 16. 

To obtain the velocities in the viscous layer 
when the wave consists of a long intumescence, 
we require a general solution with the value of u 
given by eq 15. The general solution may be 
constructed from the elementary solution, eq 16. 
Multiplying its right-hand member by F(a)da and 
integrating between the limits, + 00 and -co, we 
obtain 

U - U ^~ry^) e ' V COs[K(a-x)-2J™~\da, 

which still is a solution of eq 13 and satisfies the 
second condition of eq 14. Since the permissible 
values of k are positive numbers, we finally take 



u—u { 



cos 



1 /•«, r+ m -zJ~ 

--- dK\ F(a)e V2 
irjo J-00 

) L <oL-x)-Z^™ v ^da, 



(18) 



and this is the desired solution. For, when we 
make 2=0, u=0, we have 

I /-co p+«> 

— U = dK F(a) cos [k((x— x)]da. 

TTJO J-00 

and, therefore, according to the rule of the Fourier 
integral, it gives u =F(x), indicating that the 
first condition in eq 14 is also satisfied. 
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Although eq 18 gives the distribution of veloci- 
ties in the boundary layer for the most general 
case, it is desirable to transform the solution into 
a form more directly manageable for calculation. 
That is, it is better to change the double integral 
into a single integral by effecting an integration 
with respect to k. With a view to making use of 
the two well-known definite integrals, 



and 



f 



e gs cos (2s 2 -\-qs)sds = 0, 



(19) 



r 



T ~~~± 



e~Q s cos (2s 2 — qs)sds=^ qe \ (20) 

o 



U — U a 



- r F(a)da r 

F(a)da 

Kjx Jo 



eq 18 is broken into two parts, as follows: 

, [kco 



cos 



In the first double integral in eq 21 a^x, and this 
suggests the transformations 



and 



k(ql — x) = — 2s 2 

VKCx) 
2v =<l8 > 



I k(<x— x)— 2-y^J«?« 

VK03 _ 

cos k(ol—x)—z«Jy \dK. 
or putting q=2(3, 



(21) 



u—u a 



dn=— — ds. 
x—a 



Introducing these transformations in the first 
double integral of eq 2 1 , the integral changes to 



f* F(a)da r° 
J-oo X— a Jo 



e~y s cos (— 2s 2 — qs)sds, 



which is the desired solution and may be compared 
with the Duhamel solution of the heat-conduction 
equation [2], 

V. Loss of Energy in the Laminar Layer 

As it is assumed that the laminar layer in the 
vicinity of the solid boundaries is thin and of 
considerable length, it suffices to express the dis- 
sipation function in the form 



and therefore vanishes by virtue of eq 19. In the 
second double integral of eq 21, a^x, and this 
suggests the transformation 



*->($ 



K(a — x) = 2s 2 , 

VKOO 
2 v =q - S ' 



dn=— — ds. 
a — x ' 



where ^ is the viscosity of the liquid. Consider a 
rectangular column of indefinite height, the cross 
section being of unit length and of width dx. De- 
noting the rate of dissipation of energy in the 
column by d/dt bE\, 



d 
dt 



ol — x=- 



v q 2 



and 



u- 



(SE 1 ) = ^dx^(^J dz. (23) 

Now, since u is independent of z, 

Integrating the right-hand member by parts, we 
have, since u—u =—u when 2=0, and u—u =0 
when z= oo } 



7 2co2 2 , 

da= — — -o- dq. 
vq s * 



Introducing these transformations in the second 
double integral of eq 21 and making use of eq 20, 
the final result is 
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where (&u/dz) is the gradient of velocity at the 
wall. In view of eq 13, 



-(*)+;&* ^r^ 



/;©'- 



Substituting in eq 23, 

(SEt) = fxdydxu \^J Q + 

, . d f 00 (u—u ) 2 , 






When this expression is integrated between the 
limits x= + oo and x= — », the result gives the 
rate of energy dissipation in the entire length of 
the laminar boundary layer. Denoting this loss 
by dEi/dt, the integration gives 



dEi 

dt ' 



:+^j_ + j^(^) ( 



dx. 



(24) 



since at the infinite distances, x= «> and x = — oo , 
the velocities w and u both vanish. This formula 
is from Boussinesq [1]. The interpretation of the 
formula is that for the determination of the rate 
of dissipation of energy in the entire wave it is 
sufficient to consider the gradient of the laminar 
velocities at the solid boundary and the velocities 
in the region outside of the boundary layers. 

This result can be extended to limited waves 
with a new interpretation. Let the posterior and 
anterior feet of the wave be at the points x x and x 2 , 
respectively. Since the thickness of the viscous 
layer at the front of the wave; that is, at point x 2 , 
is negligibly small, the expression for the rate of 
dissipation of energy for the limited wave per 
unit width of channel may now be written, since 
u is independent of z, as 

where u x is the laminar velocity at the point a*. 
The interpretation of the two integrals on the 
right-hand side of the equation is as follows. The 
first integral (double) represents the rate of energy 
dissipation in the laminar boundary layer under 
the wave. The second integral represents the 
kinetic energy left behind in the channel as the 
wave travels onward. As far as the wave is 
concerned, the energy left behind is lost, so that 
also in the case of a limited wave the integral on 
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the left-hand side represents the rate of loss of 
energy of the wave per unit width of channel. 

VI. Gradual Extinction of Rectangular 
Waves 

The rectangular wave constitutes the simplest 
example of a translation wave, and the application 
of the above theories to obtain the expression for 
the gradual damping of the wave with distance 
involves no difficulty. In a rectangular wave the 
superelevation is constant, and the wave length 
is finite. 

Denote the constant superelevation of the wave 
by h u the length by L, and the particle velocities 
in the potential region by V. Assuming that the 

wave is traveling in the direction of positive x, 
we have 



K. 



U = V= to jy> to = -yJgH. 



(25) 



Suppose that at time t=0, the front of the wave 
is at the station x=L (see fig. 1). Then for the 



f 777777777T7fn7Trrrr777T7T77T7T7777777n7Tn777TT77777TT77777777T, *** 

Figure 1. Surface profile of rectangular wave 

description of the function F appearing in eq 15, 
we have 

u =F(x)=0, x^L, 

Uq= F(x) = V, O^x^L, (26) 

^ =iT( x ) = o, x^O. 

Accordingly, using eq 22, we find the boundary 
layer velocities to be 

u-V=-^f P e-?dt3, (27) 

■yj IT Jf / a> 

when x lies in the range to L, and 
tt= 4f_ r ^4f e-tdt, (28) 

when #<0. Since, in the present example, the 
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rate of dissipation of the energy of the wave per 
unit width of channel reduces to 



dE, 



-*/.**• (i?), & > 



(29) 



by the use of eq 23, it is sufficient to consider only 
eq 27. Differentiating the latter equation with 
respect to z and remembering that V is independ- 
ent of z, 



ba 



V 



4KL-X) 



bz ^>jrv(L — x)lu 

Putting z=0, the velocity gradient at the solid 
boundary is 






■\TTV(L — X)/(jJ 

Substituting in eq 24 and remembering that 
dE 1= C L VHx 



dt 



Jo ^wv(L — x)/o) 



Integrating and afterwards substituting from eq 25, 
we have finally 

^ =~ pg^H-^hlD'V >\ (30) 

at ytj- 

which is the value of the energy lost per unit time 
and per unit width of the channel. As was ex- 
plained previously, the loss is made up of two 
parts. One part represents the dissipation of 
energy in the boundary layer under the wave; the 
other part represents the kinetic energy left in the 
liquid and in the vicinity of the solid boundary 
as the wave leaves a given station. 

The energy of the wave, half potential and half 
kinetic, is given by 

E 2 =g P h 2 1 L, 

when computed for unit width of channel. The 
variation of the energy with the time is 



dE 2 



dh x 



--2gph x L-j 



(31) 



dt ~**"" 1 ~ dt 
Denoting the distance traveled by the wave by s, 

dhi_dhi ds 
dt ds dt' 



and as ds/dt represents the velocity of the wave 



Substituting this in eq 31, 



f t-MWl* f • 



(32) 



By the principle of the conservation of energy, 
the decrease of the wave energy per unit time 
must equal the dissipation of energy per unit 
time. Accordingly the right-hand members of eq 
30 and 32 are equal to each other in absolute value, 
but they are of opposite sign. With the grouping 
of the terms after the comparison is made, we 
have 



which is the differential equation giving the rate 
of decrease of the wave height with the distance 
of travel. The equation was derived by assuming 
that the channel is very wide and the effect of 
the side walls could be neglected. If the channel 
is narrow, however, the second term in the right- 
hand member must be multiplied by {1+2HIB), 
where B denotes the width of the channel. That 
is, in a channel of finite width, 



B )\H) Kg 1 '*! 



1 + 



r<i) 



H i/2 ) 



The solution of the equation is, subject to the 
assumption, h=h ol when s=0, 



H~ H e 



(33) 



where 



-£o+¥)(gru-r- <w 

Accordingly, the decrease of wave height with 
distance is given by an exponential, and a loga- 
rithmic decrement that is independent of the height 
of the wave exists. For waves traveling in water 
of constant depth, the damping decreases with 
increase of wave length, with the increase of 
channel width, and with the decrease of kinematic 
viscosity. 
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VII. Gradual Extinction of Solitary Waves and 



Denoting the superelevation of the crest by hi, 
the superelevation of the wave profile for t con- 
stant (£=0) is given approximately [3] by 



where 



h=hi sech 2 kx, 

WUi 
4ff" 



k= 



(35) 



(36) 



The dimensionless plot of the wave profile, when 
the depth of water is twice the crest elevation, 
is shown in figure 2. Using eq 9 and 12, we may 




• 1 >?>>>?>?>\??S /;?/>?>>>»;)//;//? ?))?;;j 17777777779 '777777777777777777777?. 

4 -3 -2 -I +1 4-2 +3 +4 

x/H 

Figure 2. Surface profile of solitary wave, (hi = H/2). 

write for the particle velocities in the potential 
region the expressions, assuming that the wave is 
traveling in the direction of x negative, 



u*=A sech 2 kx: A—c 



h 



(37) 



In the region of the viscous boundary layers, using 
eq22, 

u—u = 7-1 sech 2 (kx+^^j e-^dp (38) 

Differentiating with respect to z and remembering 
that u is independent of z, 

^— =H — p= sech 2 ( kx+'T ^ ) tanh ( kx+ 
ko)Z 2 \ kuz .2 jo 

W)2 V W e dp - 

The limiting value of the integral when z ap- 
proaches zero cannot be obtained directly, and a 
transformation of variables becomes necessary. 
Putting 

kx==n, 



so that 



we have 

du 
dz r 



jkoo Z 



-V: 



47 ,3 2 ^ =rfa ' 



= + 



4A 



M 



sech 2 (?i+a 2 ) tanh (n-\- 



a 2 )e a da. 
This gives for z=0 

(©r^ V?J> h2 (n+a2) tanh (ft+a2) **- 

(39) 

The rigorous determination of the definite 
integral appearing in the right-hand side is difli- 
cult, and sufficiently accurate results can be 
obtained if the indicated integration is replaced 
by the process of summation. Writing 

Y(n)= [ " sech 2 (n + a 2 ) tanh (n+a 2 )da, (40) 

the desired summation for a specific Y{n) was car- 
ried out by assigning to the interval da the value 
0.1. The results of such determinations are given 
in table 1. The graph of Y(n) is shown in figure 3. 
Since Y{n) is proportional to the velocity gradient 
at the channel bed, the curve may be interpreted 
as giving the variation of this velocity gradient 
along the wave. The velocity gradient at the 
wall reverses its sign at n=— 0.55. The maxi- 
mum value of the velocity gradient at n = 0.35. 

Table 1. — Numerical values of the function Y(n) 



11 


Y(n) 


n 


Y(n) 


4.0 


0. 0000 


-0.2 


0. 1858 


3.0 


.0066 


-.4 


.0812 


2.5 


.0167 


-.6 


-.0187 


2.0 


.0440 


-.8 


-.1017 


1.6 


.0916 


-1.2 


-.1495 


1.2 


.1725 


-1.6 


-. 1408 


0.8 


.2988 


-2.0 


-. 1132 


.6 


.3526 


-2.5 


-.0820 


.4 


.3798 


-3.0 


-.0598 


.2 


.3639 


-4.0 


-.0355 


.0 


.2975 


-5.0 


-.0243 
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Differentiating witnrespect to t, and remembering 
that 

we find forTthejrate of decrease of the energy of 
the wave 



dt ' 



'Va 



pflPhl' 



dhi 

ds 



(43) 



Figure 3. Variation of velocity gradient at wall with dis- 
tance from center of wave. 

For the case when the wave height is half of the 
undisturbed water depth, the reversal occurs at 
the point, x=—0.9H. The maximum occurs at 
the point, x=0.57H; that is, ahead of the wave 
crest. The crest position is given by x=0 or 
n=0. 

In view of the substitutions from the eq 37, 39, 
and 40, and the transformation, kx=n, the formula 
for the rate of dissipation of energy, eq 23, reduces 
to 



By the principle of conservation of energy, the 
right-hand members of the eq. 42 and 43 are equal 
in absolute value, but they are of opposite signs. 
Making the comparison and rearranging the terms, 
we obtain 



f**VV^V-± / v d ( s \ 

\HJ KB; 3\ g 1/2 H* /2 a \HJ 

This equation is derived by assuming that the 
channel is very wide, and, therefore, the effect of 
the vertical side walls can be neglected. If, 
however, the channel is narrow, the second term 
in the left-hand member must be multiplied by 
[l-\-(2H/B)], where B denotes the width of chan- 
nel. That is, in a channel of finite width, 



dE l _ 4fxA- 



dt 



-yjirk 



Vt r." nn) ^ ^ (ir<i) + K i+ -#)V^^<A)=°- 



which may be written also as 



dE } _ 
dt J. 



where 



i 4^1 2 A 



■N, 



(41) 



The integration gives, denoting the initial value 
of the crest height by h ]0 , 



sech 2 ndn. 



where the coefficient of damping has the value 



The method of summations, selecting n=0.1, 
gave for iV the value 0.316, or approximately 1/tt. 
Using the value of k from eq. 36 and A from eq. 37 
and remembering that w=^JgH, the above ex- 
pression for the rate of dissipation of energy per 
unit width of channel becomes 



K =n( 1 + ? DV?7 



£/3/2 



(45) 



dEi_ 

dt '' 



=4T- 8/2 (|Y'V / y /4 Ai /4 . (42) 



As regards the energy of the solitary wave for 
unit width, we have 

E 2 =-^g{Hh l )^. 



This gives the variation of the crest height of a 
solitary wave with the distance traversed. 

Comparing eq 33 and 44, it is seen that the ex- 
tinction of a limited wave takes place in a manner 
very different than it does with an unlimited 
wave. Although in the former case the extinc- 
tion is characterized by a logarithmic decrement 
independent of height, this is not true of the 
solitary wave. The reason for this is to be found 
in the fact that the effective wavelength of a soli- 
tary wave is a function of the crest height. With 
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the decrease of the crest height, the effective 
length increases [4]. We may now accept the 
general statement that a logarithmic decrement 
independent of the wave height exists only for a 
wave that maintains a constant wavelength during 
extinction. Thus the damped oscillations of a 
liquid in a U-tube, of a standing wave in a rec- 
tangular basin and the damping of progressive os- 
cillatory wave are characterized by a logarithmic 
decrement independent of the magnitude of the 
oscillations. 

VIII. John Scott Russell's Experiments on 
Solitary Waves 

To test the validity of the theoretical results 
above given, John Scott Russell's experimental 
data on solitary waves will be considered [5]. Al- 
though his tests were made to determine the law of 
velocity of propagation of solitary waves, the 
representation of data is such that they can be 
used also to determine the law of damping of 
waves. The experiments were made more than 
a century ago in a wooden rectangular channel, 
1 ft. wide and 20 ft. long. Altogether some 55 
tests are shown for various water depths and for 
various methods of wave generation. 

Regarding the tests, Scott Russell states: "The 
original experiments are themselves given at the 
end of this paper for the purpose of enabling 
myone who may be disposed to make use of them 
or any future purpose, either of framing or testing 
a theory, to make use of them much in the same 
way as if he had himself made the experiments." 
The farsightedness and disinterestedness of the 
investigator is evidence of a scientific attitude of 
he highest type and are cited by the present 
writer, to whom Russell's data have proven in- 
valuable, as a model that other investigators might 
well follow. 

Unfortunately the temperature of the water of 
the tests is not stated. Because of this, the de- 
sired comparison between theory and experiment 
will be made by deducing the temperatures of 
water from the observed damping of the waves 
and comparing these values with the prevailing 
air temperatures on the dates of the experiments. 

In consulting the original data, it was seen that 
the tests could be grouped together and average 
values obtained for individual water depths H. 
The method of arriving at average values may be 



explained by considering the measurement for H= 
17.9 cm. It was thought that five individual 
tests with this depth ought to be sufficient. The 
original data furnish values of h\ the height of the 
wave, and s, the distance travelled. From these, 
the ratios h l jH and s/H were formed and plotted. 
Since in these tests the initial heights A 1( j were not 
the same, the plotting gave separate curves. 
Keeping one of the curves fixed, the others were 
shifted by different amounts parallel to the axis of 
abscissas so as to give a more or less congruent 
picture. The positions of the observed points after 
the shif tings are shown in figure 4. The centers of 
gravity of the points within selected increments of 
the abscissas were obtained and are given in table 
2. Similar determinations were made for other 
depths and these are given likewise in table 2. 

Table 2. Extinction of solitary waves; average values from 
Russell's experiments 



4 

77=2.69 em 


77=7.61 cm 


7/= 10.02 cm 


s/IT 


h^IF 


s/H 


hi/If 


8/H 


hxlll 


367 


0. 420 


66 


0.380 


117 


0.354 


535 


265 


255 


.305 


424 


.286 


741 


.137 


397 


.265 


693 


.212 


962 


.112 


605 


.197 


930 


.168 


1,060 


.087 


790 


.152 


1,156 


.128 


1,170 


.057 


950 


.130 


1,398 


.099 


1,300 


.050 


1,178 


. 097 


1,597 


.076 


1,502 


.055 


1,470 


.067 


1,742 


.060 


1,740 


.035 


1,780 


.053 


1,984 


.047 


2,100 


.020 


2,270 


.043 


2,503 


.030 


Waves 27, 28, 29, 


Waves 25, 26, 40, 


Waves 1, 2, 4, 7, 


31, 33, 34, 35 


41,42 


8, 12, 13, 15, 19, 21 


H= 12.92 cm 


77=15.80 cm 


77=17.90 cm 


s/II 


hxlH 


s/H 


hxlll 


5/77 


hxlll 


114 


0.342 


84 


0.235 


96 


0.131 


388 


.262 


351 


.181 


291 


.110 


623 


211 


525 


.148 


493 


.097 


765 


.174 


667 


.128 


683 


.087 


907 


.148 


870 


.095 


889 


.071 


1,088 


.125 


1,085 


.076 


1,097 


.063 


1,325 


.098 


1,278 


.064 


1,288 


.048 


1,672 


.068 


1,435 


.054 


1,494 


.045 


1,970 


.053 


1,605 


.049 


1,706 


.041 


2,515 


.027 


1,790 


.042 


1,910 


.032 


Waves 43, 44, 46, 

47 


• 
Waves 48, 49, 50 


Waves 51, 52, 53, 
54,55 



The data for 77=17.9 cm from table 2 are 
plotted in figure 5. The smooth curve drawn 
through the points gives the adjusted values of 
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Figure 4. ' Data from Scott Russell's experiments on solitary waves. 



the wave heights. The quantities h x /H and s/H 
were read from the smooth curve. Values of the 
fourth root of H/h\ were determined and are 
plotted in figure 6 with a convenient shift in 
the origin of s/H. The remaining data from table 
2 were treated in a similar manner, and the re- 
sulting curves are shown in the figures 6 and 7. 

The distribution of the points is linear, a cir- 
cumstance that accords very well with eq 44. The 
slopes of the lines determine the coefficients of 
damping, K. The viscosities of the water can be 
obtained from eq 45, using the observed values of 
K, the corresponding values of H, and putting B, 
the width of flume, equal to 30.48 cm. The indi- 
vidual determinations of water temperature are 
shown in table 3, and the average value is found 
to be 15.8° C. 



0.15 



0.10 



0.00 



Figure 5. 
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Decrease in height of a solitary wave with dUtanrp. 
traversed. 
H= 17.9 cm. 
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s/h 
Figure 6. Coefficient of damping of solitary waves. 

O, H=2.69 cm, #=9.0x10-"; □, H=7.61 cm, #=4.85x104; 
A, H= 17.90 cm, if=3.85xl0 4 . 
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Figure 7. Coefficient of damping of solitary waves. 

O, H= 10.02 cm, #=4.67x10-4; □, H= 12.92 cm, #=4.25x10*; 
A, H=15.80 cm, #=4.25xl0<. 



Table 3. 



Determination of water temperatures in Russell's 
experiments 



H 


# 


V 


e 


Cm 






°c 


2.69 


9. 00x10-4 


0. 01168 


14.2 


7.61 


4.85 


. 01046 


18.6 


10.02 


4.67 


.01133 


15.4 


12.92 


4.25 


.01132 


15.4 


15.80 


4.25 


.01215 


12.7 


17.90 


3.85 


. 01048 


18.5 



This average temperature of 15.8° C can be 
compared with the prevailing mean air tempera- 
tures at the time of the experiments. The Russell 
experiments were made in Edinburgh during the 
period, August 3 to 8, 1837. It may be said that 
when water has been introduced from an outside 
reservoir into a wooden flume of small depth in- 
side a laboratory, during the latter months of 



summer, the working temperature of the water 
of the flume is best approximated by the morning 
air temperature, rather thai) fche minimum or 
maximum daily air temperatures. Through the 
courtesy of the Meteorological Office, Air Ministry, 
Edinburgh, the air temperatures for the dates 
corresponding to Russell's experiments were com- 
municated to the Liaison Officer, Meteorological 
Office, U. S. Weather Bureau, Washington, D. C. 
The data later were transmitted to the author by 
the Chief of the Division of Climatological and 
Hydrologic Services, U. S. Weather Bureau. The 
author wishes to express his appreciation of these 
favors. The Royal Society of Edinburgh have 
permitted a copy to be made of the pertinent part 
of Mr. Adie's original register for 1837, and the 
part of the data relating to the temperatures is 
given in table 4. It is the opinion of the Edin- 
burgh Meteorological Office that the morning and 
the evening readings of the dry bulb were taken 
at 10 a. m. and 10 p. m. Accordingly, the temper- 
atures described as the morning temperatures 
should be close to the water temperatures of the 
Russell experiments. The average of the morning 
temperature for the 7 days is about 14.6° C. This 
agrees reasonably well with the average temper- 
atures obtained from the analysis of the wave 
data; that is, 15.8° C. 

Table 4. — Edinburgh temperatures, August 1837 



Day of 
month 


Morning 


Evening 


Mini- 
mum 


Maxi- 
mum 




o F 


o F 


• o F 


° F 


1 


60 


53 


49 


57 


2 


55 


53 


50 


55 


3 


61 


55 


48 


59 


4 


52 


49 


48 


57 


5 


54 


53 


44 


58 


6 


56 


54 


46 


60 


7 


60 


56 


44 


66 



The favorable agreement between the two tem- 
peratures and the fact that the inverse of the fourth 
root of the wave height varies linearly with the 
distance of the travel of the waves are proofs 
enough to show that the present theory of viscous 
action is sufficient to evaluate the wave dissipa- 
tion for the water depths considered. The next 
step in the investigations of this type should be 
to determine what depths of water will result in a 
type of dissipation other than that caused by 
viscosity. 
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IX. Application to Model Harbor Tests 

To the extent that the theoretical formula for 
the damping of solitary waves is satisfactory, at 
least when the waves are moving in shallow water, 
we can apply the results of the above theory to 
model tests of wave action. The models are of 
the type in which the depths of water involved 
are small and in which the pattern of the individual 
waves can be approximated as a succession of 
solitary waves. No serious error is made by re- 
placing the ordinary pattern of the shallow-water 
oscillatory waves by a succession of solitary waves 
if this is done for the purpose of evaluating the 
energy dissipation. The method has the advan- 
tage that the wavelengths do not enter into the 
computations, so that the comparisons may be 
made on the basis of wave height only. 

One problem that most readily suggests itself 
and is susceptible to approximate analysis is the 
model test of wave action in a harbor. Ordinarily 
the heights of the waves approaching a model 
harbor are determined on the basis of Froudian 
scales after the respective heights of the prototype 
waves have been selected. For the actual model 
these values of the initial wave heights must be 
increased over those determined by the Froude 
law in order to correct for the energy dissipation 
in the model harbor. In the prototype the dissi- 
pation is not large, so that it can be ignored, and 
this fact facilitates the determination of the proper 
model wave heights. The exact determination of 
the desired model wave heights in the approaches 
to the model harbor depends on the length and 
depth proportions of the harbor. The evaluation 



of the energy dissipation may be made step by 
step for various consecutive segments of the model. 
The average of the time rate of dissipation for the 
entire harbor area can then be evaluated. The 
height of the approaching waves is increased to 
such an extent that the increased flow of energy 
just balances the energy dissipation in the harbor. 



The author expresses his appreciation to Herbert 
N. Eaton, Chief of the Hydraulics Section, for his 
editorial review of the paper and for constructive 
suggestions; also to George W. Patterson, his 
former colleague, for advice and the review of the 
mathematical methods used in the paper. 
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